Perturbative Renormalizability of Chiral Two Pion Exchange 
in Nucleon-Nucleon Scattering: P- and D-waves 
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We study the perturbative renormalizability of chiral two pion exchange in nucleon-nucleon 
scattering for p- and d-waves within the effective field theory approach. The one pion exchange 
potential is fully iterated at the leading order in the expansion, a choice generating a consistent 
and well-defined power counting that we explore in detail. The results show that perturbative 
chiral two pion exchange reproduces the data up to a center-of-mass momentum of fccm ~ 300 MeV 
at next-to-next-to-leading order and that the effective field theory expansion convergences up to 
fccm ~ 350 MeV. 
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I. INTRODUCTION 

The nature and derivation of the nucleon-nucleon inter- 
action is probably the central problem of nuclear physics. 
The modern point of view is that any serious theoretical 
formulation of the nuclear force should be grounded on 
quantum chromodynamics (QCD), the fundamental the- 
ory of the strong interactions. However direct calcula- 
tions on the lattice are still unavailable at the physical 
pion mass, despite promising results which are beginning 
to appear in the two [1,, ,2] and three ^ nucleon sec- 
tors at larger pion masses. A different, more indirect 
path is provided by the effective field theory (EFT) for- 
mulation of the nuclear forces in which the pion 
exchanges are constrained by the requirement of broken 
chiral symmetry, the main low energy manifestation of 
QCD. This approach exploits the well-known separation 
of scales appearing in the two-nucleon system with the 
purpose of generating a systematic and model indepen- 
dent low energy expansion of the nuclear potential and 
the scattering amplitudes. In addition, nuclear EFT may 
potentially bridge current lattice QCD calculations with 
physical results by means of chiral extrapolations. 

In the nuclear EFT proposed by Weinberg 0, [T^] the 
nucleon-nucleon potential is expanded as a power series 
(i.e. a power counting) in terms of the parameter Q/Aq 

VMr) = V^'^r) + V'^'\r) + V'^^\r) + O (^|^) , 

(1) 

where Q represents the low-energy scales of the system, 
usually the momentum p exchanged between the nucle- 
ons and the pion mass m^r , while Aq stands for the high- 
energy scales that are not explicitly taken into account in 



the theory, e.g. the rho meson mass rup. The full poten- 
tial is then iterated into the Schrodinger or Lippmann- 
Schwinger equation [lll-l33|. or more recently within the 
nuclear lattice approach j34l - l42l | 0, in order to obtain 
theoretical predictions. Regarding the notation, order 
Q*^ will also be referred to as leading order (LO), order 
as next-to-leading order (NLO), order as next-to- 
next-to-leading order (N^LO), and so on. 

Among the advantages of the Weinberg prescription we 
can count that it observes the non-perturbative nature of 
the nuclear forces and naturally fits into the traditional 
nuclear physics paradigm of employing a potential as the 
basic calculational input. In particular, the construction 
of potential descriptions of the nuclear force is conceptu- 
ally straightforward, as the procedure only involves the 
expansion of the chiral potential up to a given order and 
the determination of the free parameters of the theory by 
a fit to the available low-energy scattering data. Even- 
tually, at high enough order, in particular N'^LO, this 
procedure leads to chiral potentials that are able to repro- 
duce the two-nucleon scattering data below a labora tory 
energy of Eu^ = 300 MeV with a x^/d-o.f. ^ I 
a remarkable degree of success which justifies the popu- 
larity of the Weinberg prescription. 

However, the theoretical basis of the Weinberg count- 
ing is debatable and, as we will argue, far from robust. 
In particular we can identify two serious problems: (i) 
the counterterms included in the Weinberg counting are 
not enough as to renormalize the amp litudes (either at 
LO [ll,[2ll,[23 or NLO/N^LO ^,^) and (ii) it is not 
clear whether the scattering amplitudes obey the power 
counting of the chiral potential as a result of the full it- 
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^ It should be noted that lattice EFT implementations of the Wein- 
berg counting only iterate the order piece of the interaction, 
while the subleading pieces generally enter as perturbations. The 
C2(p^ + p''^) short range operator is an exception as it is pro- 
moted from order to order to form an "improved" LO 
interaction. 



2 



eration process. In addition, there is the issue of the 
chiral inconsistency [43| which inspired the Kaplan, Sav- 
age and Wise (KSW) counting [13, El] . Without these 
ingredients, renormahzabihty and power counting, the 
Weinberg prescription is merely reduced to a recipe for 
constructing nuclear potentials. 

The renormahzabihty problem is a consequence of the 
appearance of singular interactions in the chiral expan- 
sion of the potential: the order ly-th contribution to the 
potential behaves as y('')(r) ~ for niT^r < 1. In 

this regard the theory of non-perturbative renormahz- 
abihty developed in Refs. [2ll - [23, |46{ states that, if the 
finite range potential is singular and attractive, coun- 
terterms need to be included for removing the cut-off 
dependence. This result implies in particular the non- 
renormalizability of the Weinberg counting at any order, 
as there is always an infinite number of partial waves for 
which the chiral potential is singular and attractive. Two 
possible solutions are (i) to treat the chiral potential, in 
particular one pion exchange (OPE), perturbatively in 
partial waves with sufficiently high angular momentum, 
as advocated in Rcfs. 21, 471], or (ii) to realize that the 
infinite number of counterterms is partly an artifact of 
the partial wave projection that can be avoided by corre- 
lating the short-range operators of different channels [1^ . 
In the present work we will follow the first solution, which 
is the simplest and most natural one from the EFT per- 
spective. 

On the contrary, a repulsive singular interaction is 
completely insensitive to the value of the counterterms. 
This condition leads to unexpected consequences, in par- 
ticular the deuteron NLO disaster [23*1: at NLO the chiral 
potential turns out to be repulsive in the triplet chan- 
nel leading to the disappearance of the deuteron as the 
cut-off is removed, in clear disagreement with the ex- 
perimental evidence. In addition, the eventual change 
of sign of the chiral potential at certain orders prevents 
the formulation of a sensible power counting in a fully 
non-perturbative set-u p, a s has been recently stressed by 
Entem and Machleidt [48| . 

However, the difficulties associated with non- 
perturbative renormahzabihty are not the only reason to 
avoid the iteration of chiral two pion exchange (TPE) . A 
different, but also potentially serious issue is the possi- 
ble mismatch between the power counting in the poten- 
tial and the physical observables. As is well known, the 
Weinberg prescription defines power counting in terms of 
the NN potential, which is not an observable. The itera- 
tion of the potential can alter the relative importance of 
each component of the interaction, specially if the cut-off 
is not soft enough [3l|, [s^l , and thus destroy the 

original power counting at the level of observables. In 
this respect the preliminary results of Ref. [sj] for the 
singlet channel indicate that this may be happening for 
the typical values of the cut-off employed in the Wein- 
berg scherne, calling for a perturbative reanalysis of the 
N^LO [ll[ll,[ll and N^LO ll7j,|2^ results. 

As argued in the previous paragraphs, the combined 



requirement of renormalizability and power counting im- 
poses stringent constraints on the set of theoretically 
acceptable EFT formulations of the nuclear force. At 
leading order the non-perturbative features of the two- 
nucleon interaction requires the iteration of certain pieces 
of the interaction, minimally a contact operator in the 
^So and channels [i3,lii,IIl, and optimally the OPE 
potential in s- and p- waves (and eventually d- waves), 
plus the extra counterterms required to renormalize the 
scattering amplitude [2l|. In order to avoid counterterm 
proliferation, we should be conservative with respect to 
what pieces of the chiral potential we iterate. Finally, to 
guarantee the existence of a power counting and prevent 
the renormalization issues we have mentioned, the sub- 
leading pieces of the interaction should be perturbative 
corrections over the LO results. The modified Weinberg 
scheme of Ref. 21] fulfills these conditions and is in ad- 
dition a phcnomenologically promising approach, as can 
be inferred from its lowest order results. However, as 
we are dealing with a posteriori power counting, the pre- 
vious is not necessarily the only possible realization of 
a consistent and phcnomenologically acceptable nuclear 
EFT, and there may be other frameworks or variations 
(see Refs. [S^ [13] for two examples) that may eventually 
work as well. 

The present work explores the perturbative renormal- 
izability of chiral TPE at NLO and N^LO assuming the 
non-perturbative treatment of OPE at LO. This for- 
mulation is to be identified with the proposal of Nogga, 
Timmermans and van Kolck [2l| for constructing a con- 
sistent power counting for nuclear EFT. Regarding the 
formal aspects of the theory, we elaborate upon the per- 
turbative renormalization framework of Ref. ^55| and de- 
termine the general conditions that ensure the renormal- 
izability of the subleading corrections to the LO scatter- 
ing amplitudes. That is, we deduce the power counting 
of the contact operators. As expected the scaling of the 
counterterms is very similar to the one previously ob- 
tained by Birse 147] . As an application of the formalism, 
we compute the p- and d-wave phase shifts resulting in 
a phcnomenologically acceptable description which im- 
proves over the WeinlDerg counting at the same order [l^ . 

The general approach we follow and the perturbative 
techniques we employ, which are based on distorted wave 
Born approximation (DWBA), are equivalent to the mo- 
mentum space formulation of Refs. (sgUst} . In particular, 
the recent exploration of the '^Pq partial wave by Long 
and Yang [s^] leads to conclusions that are equivalent to 
the ones obtained in the present work for the aforesaid 
wave. The role of perturbative chiral TPE has also been 
studied in the "deconstruction" approach of Refs. [ssl - 
[gq I] , which analyzes the scaling of the short range inter- 
action for most of the uncoupled s-, p- and d- waves within 
the power counting derived from Ref. [2l[. In addition 
deconstruction provides interesting insights on important 
aspects of the theory such as the determination of the ex- 
pansion parameter of the EFT, which we will employ in 
the present work. 
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The article is structured as follows: in Sect. [IT] we 
review the modifications induced by the requirement 
of renormalizability to the Weinberg power counting at 
leading [^H and subleading (5^ orders. In Sect. lIIII we ex- 
plore in detail the divergences related to the perturbative 
treatment of chiral TPE within the DWBA to determine 
the power counting of the subleading contact operators 
at NLO and N^LO. The p- and d-waves phase shifts re- 
sulting from this renormalization scheme are presented in 
Sect. HVl Finally, we discuss the results and present our 
conclusion in Sect.|Vl In the appendix [XI we explain the 
technical details behind the LO non-perturbative renor- 
malization. 



II. MODIFICATIONS TO THE WEINBERG 
POWER COUNTING 

In this section we want to explain in an informal man- 
ner why renormalizability requires certain modifications 
to the Weinberg counting. The starting point is the 
work of Nogga, Timmermans and van Kolck [211] , which 
pointed out, by means of a thorough numerical explo- 
ration, that cut-off independence of the leading order 
nucleon-nucleon scattering amplitude requires new coun- 
terterms in partial waves for which the LO tensor force 
is attractive. This precise pattern is naturally explained 
within the non-perturbative renormalization framework 
developed in Refs. [22| - [23 |. which proved that the inclu- 
sion of a specific number of counterterms is a necessary 
and sufficient condition for the renormalizability of at- 
tractive singular potentials. Curiously, the failure of the 
Weinberg scheme is not accidental: the existence of sin- 
gular interactions in the EFT expansion of the chiral po- 
tential is just a consequence of power counting itself. We 
will illustrate this point in Sect. Ill Al For exemplifying 
the ideas behind the non-perturbative renormalization of 
singular potentials and the LO modifications of the power 
counting 21-^], we will consider the particular case of 
the '^Pq partial wave in Sect. Ill Bl 

Nevertheless, the purpose of the present work is to ex- 
plore the extension of the Nogga et al. J21|] counting to 
subleading orders. As already conjectured in Ref. |2l| . 
the most adequate approach for formulating a canoni- 
cal EFT expansion of the amplitudes is the perturbative 
treatment of the higher order corrections to the chiral 
nuclear potential and the addition of a certain number 
of counterterms to guarantee the renormalizability of the 
approach. The power counting for the subleading contact 
interactions, that is, the order of the promoted coun- 
terterms, was constructed for the first time by Birse in 
Ref. (see also Ref. 61 1 for a more accessible presen- 
tation). More recently, Ref. [11] proved the viability of 
the previous proposal and found some deviations from 
the original power counting of Ref. 47[. Incidentally, 



the subleading power counting deviates from the original 
Weinberg prescription even in partial waves for which the 
leading order counting remained unaltered in the Nogga 



et al. proposal [2l|. The best example is provided by 
the ^5*0 partial wave: although at LO both Weinberg 
and Nogga et al. predict the same number of countert- 
erms, at NLO and N^LO the counterterms predicted by 
Weinberg are not enough as to renormalize the scattering 
amplitude in this partial wave. Based on the related dis- 
cussion of Ref. j5l|, we will explain the power counting 
of the singlet channel in Sect lll Cl 



A. Singular Potentials 

The appearance of singular interactions in the chi- 
ral potential is the fundamental reason which requires 
the power counting modifications suggested by Nogga et 
al. |2l| . Curiously, singular potentials are a consequence 
of the existence of a power counting in the potential it- 
self. The argument for the previous statement is rela- 
tively straightforward. We first consider the scaling of 
the order v-th contribution to the chiral potential, which 
according to power counting is given by 



(2) 



where Q represents the light scales of the system, like 
the pion mass or the relative momentum of the nucleons, 
while Aq stands for the heavy scales, for example the 
mass of the p meson. 

It is interesting to notice that the two aforementioned 
light scales play a very different role: while the pion mass 
m^r is fixed, the momentum of the nucleons can vary, 
and in particular it can take values that are much larger 
than the pion mass, that is, \q\ ^ m^. In this regime 
the momentum space representation of V^'^'^ fulfills the 
approximate scaling law 



(3) 



which is just a restatement of Eq. After Fourier 

transforming, this translates into the following scaling 
relation in coordinate space 



(4) 



which is valid for m7r|?^ <C 1. 

There are two kinds of solution to the scaling rela- 
tions induced by the power counting of the potential: (i) 
zero range and (ii) finite range solutions. The zero (or 
contact) range solutions are constructed from the scaling 
properties of the three dimensional Dirac delta, that is 



(5) 



which corresponds to a solution of the scaling equation 
for 1/ = 0, yielding 



(6) 



where the subscript c indicates the contact range nature 
of the contribution. Higher order contact range solutions 
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can be easily obtained by adding derivatives to the Dirac 
delta, although due to parity conservation only an even 
number of derivatives can appear. That is, new contact 
terms enter only at at even orders, v ~ 2n, as is well 
known. 

The finite range solutions of the scaling relations are 
fairly straightforward and take the form 



1 



(7) 



where the subscript p is used to indicate the finite range. 
This form corresponds to the observed degree of diver- 
gence of the pion contributions to the chiral potentials 0- 
As we will see, the renormalizability problems arise from 
the non-perturbative interplay between contact and fi- 
nite range operators, which modifies their scaling prop- 
erties at the level of observables, thus changing the power 
counting of the operators. 



B. Leading Order Modifications of the Counting 

Weinberg power counting implicitly assumes that the 
counterterms appearing at each order in the chiral ex- 
pansion are capable of renormalizing the scattering am- 
plitude. However, this is not the case, as has been repeat- 
edly discussed in the literature [2l|, [H, [i^ . In par- 
ticular, Nogga et al. [2lj found numerically that even at 
leading order the counterterms of the Weinberg counting 
do not render the phase shifts unique: several countert- 
erms should be added in the ^Pq, ^Pi — ^Fi channels and 
eventually the channel, depending on the cut-off re- 
gion under consideration (see also the related comments 
of Ref. [11)0. 

The ^Pq channel probably provides the best example to 
exemplify the ideas behind the Nogga et al. proposal [2l| . 
We start by consider this partial wave in the Weinberg 
counting, where we can explicitly check the appearance 
of a strong cut-off dependence of the ^Pq phase shift at 
LO . In the Weinberg scheme, the wave function in this 
channel can be described by the following Schrodinger 
equation 



,(0)" 



where Vsp' (r) is the order chiral potential. As the po- 
tential diverges as Xjr^ at short distances, we include a 
regularization scale which serves as a separation scale 
between the unknown short range physics and the known 



uu(r) = euf{r), (8) 



The inclusion of light but static degrees of freedom, like the A 
isobar in the small scale expansion 62], can alter however the 
previous scaling. 

Notice however that the recent application of the N /D method to 
nuclear EFT [ssl does not require the inclusion of an additional 
counterterm in the ^Po or the ^D-i partial waves. 



long range physics. In particular, we only consider the 
chiral potential to be valid for r > r^. For radii below 
the regularization scale, r < r^, we do not know how 
the system can be described, but it is reasonable to as- 
sume that at short enough distances the behaviour of 
the wave function will be dominated by the centrifugal 
barrier, that is. 



(9) 



for r < Tc- The previous assumption can be effectively 
translated into a logarithmic boundary condition for the 
Schrodinger equation at r — rc 



(10) 



corresponding to the radial regulator employed in 
Ref. [2J|. By integrating the Schrodinger equation from 
r = Tc to r — ^ cx) with the previous initial integration 
condition, we can extract the phase shifts by matching 
to the asymptotic form of the wave function at large dis- 
tances, which is given by 



cot ji (kr) 



yi{kr) 



(11) 



where (Jgp^ is the LO ^Pq phase shift, and ji(a::) = xji{x), 
yi{x) — xyi{x), with ji{x) and yi(x) the Spherical 
Bessel functions. 

The results for the "^Pq phase shifts in the previous reg- 
ularization scheme are depicted in Fig. ([1]). The phase 
shifts, which have been computed for a series of cut-off 
radii ranging from = 0.6 fm to 1.6 fm, show a remark- 
able cut-off dependence not anticipated by the Weinberg 
counting. The formal reason for this cut-off dependence 
lies in the details of the short range {mT^r ^ 1) solution 
to the Schrodinger equation Eq. ([8]), which reads 



u^^\r) cx (fcaa) 



^3, 



3/4 






j sin 


\ — + V 
V r 



{12) 



where is a length scale related to the strength of the 
tensor force in this channel. As can be seen, the determi- 
nation of the wave function near the origin requires the 
determination of the phase ip. The boundary condition, 
Eq. (fTO|l . just chooses a different, non-converging value 
of Lf for each r^, and therefore does not yield unique re- 
sults. In particular, there is not a well defined limit of 

for Tc — 0: the value of i/j simply oscillates faster and 
faster on the way to the origin. 

From the EFT point of view, the previous cut-off 
dependence means that the a priori estimation of the 
strength of the contact term for the ^Pq channel was not 
correct. Otherwise, the cut-off dependence would have 
not appeared until radii of the order of the breakdown 
scale of the theory AqTc ^ 1. However, as can be seen in 
the left panel of Fig. ([1]), the cut-off dependence manifest 
for scales of the order of to^Tc ^ 1. This indicates the 
necessity of including a new short range operator in this 
channel at LO. 
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FIG. 1. (Color online) Leading order phase shifts for the '^Po channel in the Weinberg counting (left panel) and in the modified 
Weinberg scheme of Nogga, Timmermans and van Kolck 21] (right panel). As can be seen, the ^Po phases develop a strong 
cut-off dependence in the Weinberg counting, which can be eliminated by the promotion of a counterterm as advocated in 
Ref. dH. 



By including the new counterterm, the previous uncon- 
trolled cut-ofF dependence of the '^Pq phase shift disap- 
pears, as can be appreciated in the right panel of Fig. (111 . 
We have employed the regularization scheme of Refs. |22| - 
[23 |. in which to unambiguously determine the LO phase 
shifts in the ^Pq partial wave we fix the value of the scat- 
tering length 0, see the appendix |X] for the details. Of 
course, there is still a residual cut-off dependence that is 
however mostly harmless owing to the existence of a clear 
convergence pattern in the rj, — > limit. That is, we can 
interpret the residual cut-off dependence as a higher or- 
der effect. For the regularization employed in this case, 
the cut-ofF dependence of the phase shifts can be explic- 
itly computed by employing the methods of Ref. |64| . 
yielding the result 



,(0) 



{rc) 



k 



1/2 
3 ■ 



(13) 



which serves as a formal confirmation of the existence of 
the Tc — 7- limit. 



C. Subleading Order Modifications of the Counting 

The modifications of the power counting induced by 
renormalizability are not confined to leading order. Fur- 
ther deviations occur at subleading orders, which can 
be identified by analyzing the cut-off dependence of the 



* Of course, using the term scattering length for a p-wave is lan- 
guage abuse. The proper name is scattering volume, as it has 
dimensions of [length] ^ . A similar comment would apply in the d- 
wave case, in which the so-called scattering length has dimension 
of [length]^ and henceforth is a (5-dimensional) hypervolume. 
However, for notational simplicity we will always use the term 
scattering length for the coefficient related to the low energy be- 
haviour of the i-wave phase shift, <5i(fe) = —a; A:^'+^ -|- C'(fc^'+^). 



phase shifts, even though these orders are included as 
perturbations. In this case, the most direct example is 
the ^S'o channel, which follows the Weinberg counting at 
LO, but diverts from it when subleading corrections are 
included. At orders and the Weinberg counting 
dictates a total of two counterterms for the ^Sq channel, 
which in momentum space take the form 

{p\vt'^%')=C^^'^+C^^\p'+p'^). (14) 

In terms of renormalization, this counterterm structure 
is equivalent to fixing two observables, for example the 
scattering length and the effective range of the ^ 5o chan- 
nel. For making such a calculation we will follow the 
formalism of Ref. 55| , but adapting the number of coun- 



terterms to two. In such a case, we obtain the phase 
shifts of Fig. dl]). 

As can be appreciated, there is a strong cut-off depen- 
dence in the results of Fig. The reason can be found 
by analyzing the cut-off behaviour of the perturbative 
corrections to the phase shifts. According to Ref. [s^, 
the Q'^ contribution to the phase shift is proportional to 
the integral 



(5('')(fc;rc) cx 



A* 



drV^''\r)u^^^ (r 



(15) 



where u'^'^ {r) is the LO wave function in the Wein- 
berg counting, which can be described by adapting the 
Schrodinger equation for the ^Pq channel, Eq. ([TU]), to 
the ^5*0 case, and V^"'' is the order Q" contribution to 
the chiral potential. It should be noted that 5^"^ scales 
as Q'^'^^ : even though nominally it is of order Q'^ . The 
degree of divergence of 5'^^^ can be easily determined by 
considering the Taylor expansion of the LO wave function 
fH, that is 



Ufc°^(r) = uo + k'^U2 + k^Ui + 



(16) 



where the behaviour of the different terms at short 
enough distances is given by mq ~ 1, M2 ^ t"^, M4 ~ r*. 
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FIG. 2. (Color online) Next-to-next-to leading order phase shifts for the ^So channel in the perturbative Weinberg counting 
(left panel) and in the Nogga, Timmermans and van Kolck scheme [IJ (right panel). By perturbative Weinberg we mean that 
the subleading pieces of the potential (that is, chiral two pion exchange) are treated as perturbations. We can appreciated 
that the cut-off dependence appearing in the perturbative Weinberg set-up can be cured by the inclusion of an additional 
counterterm. 



and so on. With the previous information, and the ad- 
ditional fact that V^"'^ ~ l/r^+'', we find that for the 
particular case v = 3 the integral in Eq. ([15]) diverges as 

/"Oo 2 

&=^^k;r^)^ drW'')(r)4°) (r) 




with a the length scale governing the divergences and 
where cq, C2 and C4 are dimensionless numbers that are 
expected to be of order one, while the terms in the dots 
are already finite. In Fig. ^ we have only fixed the fc" 
and fc^ behaviour for reproducing the low energy phase 
shifts, that is, the scattering length and the effective 
range. This means that there is still a contribution that 
diverges as k*/rc explaining the failure of the Weinberg 
counting in the perturbative context. 

Again, the solution is to modify the counting in such 
a way as to include a new counterterm. This means that 
the correct and contact interaction should be 

+ Cf)(/+p'4), (18) 



^ We have ignored the C^'^'p^p'^ operator, as it is redundant |65| . 
However, the equivalence of this operator with C4(p'* +p'^) only 
shows up either in dimensional regularization or once the floating 
cut-off is removed. For finite cut-off calculations, the difference 
between the C4 and C"^ operators is presumably a higher or- 
der effect. On the practical side the inclusion of this operator 
will surely be advantageous in momentum space treatments as 
it may help to reduce the cut-off dependence of the scattering 
observables. 



instead of the form dictated by Weinberg dimensional 
counting in Eq. The additional counterterm is able 

to absorb the fc^/rc divergence in the N^LO phase shift. 
In such a case, the new cut-off dependence is the one 
given by Fig. ([2]), in which there is still a residual cut-off 
dependence of 0{k^rc). However, the same comments 
apply as in the '^Pq channel at LO : the residual cut- 
off dependence is a higher order effect and can be safely 
ignored. 

Of course, the power counting modification of Eq. 
is well-known to arise in two-body systems with a large 
value of the scattering length [i^, IH, [H, [63| , where the 
renormalization group analysis of the C4 operator indi- 
cates that its actual scaling is rather than the naive 
dimensional expectation Q'^. In this regard, it should 
be stressed that the ^Sq phase shifts cannot be properly 
renormalized at NLO/N^LO without the inclusion of the 
C4 operator. The previous is not merely a statement 
about the existence of the Tc — >■ limit in the calculation, 
but also about the accuracy of the theoretical description 
of the system. The renormalized calculation is accurate 
up to Q^iQ^) at NLO (N^LO). On the contrary, ff we 
do not to include the C4 operator, the error of the calcu- 
lation will be 0{Q^), much bigger than expected. This 
may prove to be a si gnif icant drawback for the nuclear 
EFT lattice approach |36l - l42| , which does not include the 
aforementioned operator in the ^5*0 channel. 

Finally, the observation that the contact interaction 
of Eq. (jlSp renormalizes the scattering amplitude can be 
employed to extract the breakdown scale of the EFT in 
the singlet channel, as has been shown by the deconstruc- 
tion approach of Refs. [H, [11]. The idea is to employ 
the phenomenological phase shifts as the input of a cal- 
culation for determining the form of the contact range 
interaction up to an error of 0{Q'^^^). This represents 
in fact the complementary process of fitting the Cj^^ pa- 
rameters to the phase shifts. By considering the on-shell 
form of the contact range potential, V^'"\k) = (fcll^^^lfc). 
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Ref. [13 is able to determine that (i) the form of the 
short range interaction compUes with Eq. ([T5| and (ii) 
the hard scale governing Vc is Ao,s — 270 MeV. This fig- 
ure is much lower than the expected Aq ^ 0.5 GeV and 
translates into a slow convergence rate for the singlet 
channel, Q/Aq ~ TOtt/Acs — 0.5, emphasizing the im- 
portance of the proper renormalization of the scattering 
amplitude to minimize the theoretical uncertainty. The 
breakdown scale can in fact be appreciated in the right 
panel of Fig. where the cut-off dependence of the 
renormalized results starts to become sizable at fc ~ Aq^s, 
in agreement with the estimates of Ref. 'b^ 0. 



III. PERTURBATIVE RENORMALIZABILITY 
OF SINGULAR POTENTIALS AND POWER 
COUNTING 

In this section we study the power counting that arises 
when one pion exchange is iterated to all orders and chi- 
ral two-pion exchange is treated as a perturbation. We 
follow a more formal and detailed approach than in the 
previous section. We start by defining the power counting 
conventions in Sect. IIII Al We derive the power count- 
ing by requiring the perturbative corrections to remain 
finite in the Tc — limit, first for the uncoupled chan- 
nels in Sect. IIII Bl and then we extend the results to the 
coupled channels in Sect. IIII CI We informally comment 
on the convergence rate and the expansion parameter of 
the resulting EFT in Sect. IIII Dl Finally, we discuss the 
perturbative power counting in Sect. IIII El The results 
of this section are briefiy summarized in Table HI where 
we can check the number of counterterms required in the 
lower partial waves from LO up to N'''LO. 



A. Conventions 

We assume that the chiral potential is expanded in 
terms of a power counting as follows 




where Q (Ao) is the soft (hard) scale of the system. The 
expansion starts at order vq > —1, and we only consider 



A similar comment applies for the phase shifts computed in the 
Weinberg counting: if we inspect the left panels of Figs. JTJ 
and l(2]l, the cut-off dependence is already conspicuous at fc ~ 
100 — 200 MeV, suggesting a relatively soft breakdown scale for 
both the ^ So and ^Pq partial waves. In contrast, for the ^Pq 
channel in the Nogga et al. counting [Slj the cut-off dependence 
of the right panel of Fig. JTJ does not provide a clear indication 
of the breakdown scale. In this case the determination of the 
range of applicability of the theory is better done by analyzing 
the size the subleading order corrections to the phase shift. 



the corrections to the chiral potential up to a certain or- 
der z/max (= 2,3 in the present paper). The T-matrix 
follows the same power counting expansion as the poten- 
tial, that is 

T = T^") + O [(^)''"''^+^~) . (20) 

The scattering equations can be obtained by reexpanding 
the Lippmann-Schwinger equation 

T + VGoT , (21) 

in terms of the power counting of the potential, taking 
into account that the resolvent operator Gq is formally 
of order Q. In this convention, only contributions to 
the chiral potential which are of order Q^^ should be 
iterated, yielding the equation 

T(-i)^\/(-i) + y(-i)GoT(-i). (22) 

The corresponding dynamical equation of the higher or- 
der contributions to the T-matrix are increasingly more 
involved. 

In principle, in the Weinberg scheme the power count- 
ing expansion starts at order Q° 

V^N ^V^''^ +V^^^ +V^^^ +0{Q^). (23) 

Within the power counting conventions we follow, the 
previous means that all pieces of the potential should 
be treated as perturbations. However, it is clear that 
needs to be iterated in certain partial waves. This 
requires a promotion from order Q'^ to Q^^, that is 

1/(0) ^ , (24) 

at least in the partial waves in which either the LO 
contact operator or OPE are non-perturbative. From a 
purely quantum-mechanical point of view the iteration is 
necessary if the potential is not weak. In terms of power 
counting, a large coupling constant can be accounted for 
by identifying it with Aq/Q, which is a large number. Al- 
ternatively, we can try to explicitly identify the low en- 
ergy scale which requires the LO potential to be iterated. 
For the LO contact operator the additional low energy 
scale is well known to be the inverse of the large scatter- 
ing length of the ^Sq and ^Si — ^Di channels, as has been 
extensively discussed in the literature [3, |45|, [H, [H, [gI] . 
For justifying the iteration of the OPE potential Birse 
has proposed [i^l the identification of 

AoPE = Ann = =^ 300 MeV (25) 

as the new low energy scale. However, this choice does 
not explain the fact that OPE remains perturbative in 
the singlet channels even if iterated . A different pos- 
sibility may be to treat the inverse of the length scale 
governing the strength of the tensor force as the miss- 
ing low energy scale, giving a Aqfe that is in general a 
fraction of Ann- 
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Partial wave 


LO 


NLO 


N^LO 


N^LO 


'So 


1 


3 


3 


4 


■'Si - ^Di 


1 


6 


6 


6 


'Pi 





1 


1 


2 


'Po 


1 


2 


2 


2 


'Pi 





1 


1 


2 


■'P2 - ''F2 


1 


6 


6 


6 













1 




1 


2 


2 


2 


■'D3 - 'Gs 











1 


All 


5 


21 


21 


27 



TABLE I. Total number of counterterms per partial wave at 
LO iQ-'), NLO (Q^), N^LO {Q') and N^LO (Q") assuming 
(i) the iteration of the OPE potential at LO in the ^Si — 
^Di, 'Po, ^P2 — ^F2 and 'D2 partial waves and (ii) energy- 
dependent counterterms. The N"^LO results assume that an 
energy-independent representation of the finite-range chiral 
potential is being used. If an energy-dependent potential is 
employed, the number of free parameters will grow from 27 
to 35 (see the discussion in Sect llllEl for details). 



After the promotion of LO potential from order 
to order Q"^, we encounter a remarkable simplification 
in the dynamical equations describing the and 
contributions to the T-matrix, namely 

Tt"^) = (1 + T^-^^Gq)V^''\GoT^-^^ + 1) , (26) 

where, due to the absence of contributions to the finite 
range chiral potential from order to order , second 
order perturbation theory is not needed up to order . 
That is, the equation above is vaUd for v < b. 



B. Uncoupled Channels 

We consider in the first place the power counting ex- 
pansion of the Z-wave phase shift 



a normalization factor and the perturbative integral. 
The perturbative integral /^^'^^ is defined as 



/(^)(fc;r,) 



drV^^^r) W 



(29) 



where 



to be asymptotically (r 



is the LO wave function, which we assume 
00) normalized to 



;^'(r) — ^ cotJ^ 



yiikr), (30) 



where j;(a;) = xji{x), yi{x) = xyi{x), with j/(a;) and yi{x) 
the spherical Bessel functions. For the purpose of analyz- 
ing the renormalization of the perturbative phase shifts, 
we will assume that the normalization factor A[ (fc) 
is defined in such a way that the LO wave function is 
energy-independent at r — >■ 0, that is 



(31) 



a prescription that will make relatively easy the identi- 
fication of the necessary number of counterterms. How- 
ever, in practical calculations we will simply take 



dk 



= 0, 



(32) 



which is easier to implement than Eq. (j3ip . 

The renormalizability of the perturbative phase shifts 
can be easily analyzed in terms of the short range be- 
haviour of the LO reduced wave function and the 
contribution to the chiral potential, which determine 
whether the perturbative integral is divergent or not. 
In particular, for mr <^ 1, we have the following 



(r) 



1 



(33) 



where s describes the power-law behaviour of the wave 
function at short distances. For the perturbative integral 
we have 



^+1 



). (27) 



(34) 



The LO phase shift is calculated non-perturbatively by 
solving the reduced Schrodinger equation with the LO 
chiral potential, see Appendix [X] for details. On the 
other hand, we compute the order ly-th contribution to 

the phase shift, Sl'^'\ in the distorted wave Born approx- 
imation (DWBA), that is 

Sl''\k;r,) 2^l 



sm 



A 



/r(fc;r,) (28) 



where is the leading order phase shift, fi the reduced 
mass of the two nucleon system, k the center-of-mass mo- 



mentum, / the angular momentum of the system, A 



(-1) 



which diverges if 3 -I- — 2s > 1. The necessary number 
of counterterms can be determined by considering the 
Taylor expansion of the perturbative integral in terms of 
k^ 



(35) 



ri=0 



in which each new term is less singular than the previous 
one as a consequence of the related expansion for the LO 
reduced wave function 



Tt = 



{r)k'^ 



(36) 
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where 1^2"^ '-^ ^s+nt j^j. qt^qH radii, with t > 2. We 
can see now that the reason for choosing the energy- 
independent normahzation of Eq. pip is to have addi- 
tional power law suppression at short distances in the 
expansion of the wave function. In terms of the expan- 
sion of the perturbative integral, we obtain 

4:!.(^c)-£ ,3^/_L„. , (37) 

which converges for 3 + 1^ — 28 — 711 < 1. Eventu- 
ally, as n increases, we will have well-defined integrals. 
In this regard, we can define the number of subtrac- 
tion / counterterms as the smallest value of Uc for which 

3 + iy-2s-nct <1. (38) 

It should be noted that Uc refers to the total number 
of subtraction/counterterms required at a given order. 
For example, if ric = 1 the counterterm only affects the 
scattering length of the i^-th order of the phase shift. If 
the scattering length was already fixed at a lower order, 
this just indicates that a perturbative correction must be 
added to the lower order counterterm. 

To summarize, we have that the perturbative integral 
can be divided into a divergent and a regular piece 

//^)(fc;r,)=/[';](fc;r,)+/(';](fc;r,), (39) 

where the divergent piece can be expanded in powers of 
the squared momentum fc^, 

toik^rc) =J2lLlir c)k'\ (40) 

71=0 

with Uc the number of counterterms that renders the 
perturbative phase shifts finite. The perturbative inte- 
gral can be regularized in any specific way which we find 
convenient. A particularly simple and straightforward 
method is to add Uc free parameters to the original inte- 
gral as follows 

ric — l 

//"^^IfcirJ^ 5] Ag„fc^" + /r^(fc;re), (41) 

where the parameters x['^2n *° fitted to the experi- 
mental phase shifts. The previous parameters, which we 
will informally call A-parameters, can be related to the 
more standard counterterms if we consider a concrete 
representation of the short range physics, for example 

^S'(r-c) = ^ E CiZirc)k'- Sir - , 

(42) 

where fi{rc) = ^'^''^i^^'' is a factor which ensures that the 
matrix elements of the previous potential in momentum 
space behaves as 

{p\Vc,^\p')^p'/Y.^2nArc)k'\ (43) 

n 



for Tc — > 0. If we have chosen to work in the practical 
normalization of Eq. (|32p . the following relationship is 

obtained between the fitting parameters A2^\ and the 
usual counterterms 

For other representations of the short range physics, the 
relationship will take a more complex form. 

We can distinguish four cases: (i) the irregular solution 
of a non-singular potential (^Sq), (ii) the regular solution 
of a non-singular potential (^Pi, ^£'2), (iii) the general 
solution of an attractive singular potential {"^Pq, ^02) 
and (iv) the regular solution of a repulsive singular po- 
tential i^Pi)- In the following paragraphs we will discuss 
each of these cases in detail. 



1. Non-Smgular Potential 

Non-singular (or regular) potentials are potentials for 
which the regularity condition is fulfilled 

\im r^V{r) = 0. (45) 

r— ^0 

The regularity condition implies in turn that the short 
range behaviour of the wave function is determined by 
the centrifugal barrier, which overcomes the potential at 
short distances. For r — 0, we can define a regular and 
irregular solutions which behave as 

4;.Veg(0 - r-'+l , (46) 

«wl(0-^. (47) 

If there is no short range physics at LO, the regular so- 
lution is chosen. On the contrary, if a contact operator 
is included at LO, the LO reduced wave function will be 
a superposition of the regular and irregular solutions. 

The regular solution gives rise to the Weinberg power 
counting unaltered. If we consider the k'^ expansion of 
the regular solution, Eq. ([36)) . we have 

- . (48) 

that is, s = I + 1 and t = 2. The convergence of the 
perturbative integral requires in this case 

ly <2l + 2nc, (49) 

where stands for the number of counterterms in a given 
partial wave. Therefore new counterterms appear as ex- 
pected in naive dimensional analysis, that is 

C2n,l ^ Q''+'" . (50) 

For the particular case of the ^Pi channel, the first coun- 
terterm (Cq/Pj) is required at Q^, a second (C2.ipj) at 
Q^, and so on. 
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In contrast, the occurrence of the irregular solution at 
LO will change the power counting. The fc^ expansion of 
the wave function yields in this case 

-1". -l - , (51) 

which implies s = —I and t = 2. The finiteness condition 
reads 

2nc>v + 2l + 2, (52) 

which implies a significant deviation from naive dimen- 
sional analysis. In fact the scaling of the counterterms is 
given by 

For the ^5'o channel, the previous requires that the first 
perturbative correction to the Cg 15^ counterterm enters 
at order Q"^, while the C2„.iSo counterterms (with n > 
0) are required at order Q^"~'^. It should be noticed 
that the order assignment to the first perturbative 
correction to the Co counterterm is merely formal: this 
correction actually enters at the same order as the first 
perturbative correction to the finite range piece of the 
potential, which happens at order Q^. In fact, the 
scaling is rather an indication of the fine-tuning required 
for having an unnatural scattering length. In terms of 
Birse's renormalization group analysis (RGA) [ggI ). the 
Q^^ eigenvalue is a refiection of the unstable nature of 
the infrared fixed point associated with two-body systems 
with a large scattering length. 

It is interesting to notice that the extension of the 
power counting for irregular solutions to p- waves leads to 
the same conclusions as p-wave Halo EFT [tO], namely, 
that the first two contact operators should be iterated. 
The previous arguments imply that for a p-wave, the first 
counterterm will be of order Q^*, the second Q~^, and 
the third Q'^. This has two interpretations; first, that 
the unstable infrared fixed point is even more fine-tuned 
than in s-waves, and second, that it may be necessary to 
iterate two counterterms instead of only one. However, 
for the particular case of nucleon-nucleon scattering, we 
do not encounter this situation. 



2. Singular Potential 

The appearance of power-law singular interactions at 
LO entails substantial changes to the power counting of 
the short range operators. In particular we are interested 
in the effect of tensor OPE, which at short distances (r — >■ 
0) behaves as 

2fiVi-ll{r) ^ ±^ , (54) 

where /i is the reduced mass of the two-nucleon system 
and 03 a length scale that governs the strength of the ten- 
sor force. We have chosen a convention in which > 
and where the attractive or repulsive character of the 



potential is indicated by the explicit ± sign. The renor- 
malization of singular potentials implies that attractive 
singular interactions require the inclusion of a countert- 
erm, while repulsive singular interactions do not f22l-f2^ , 
a feature which will be reviewed in the following para- 
graph. 

For a tensor force, the terms of the expansion of the 
wave function behave for r — as 

„M)(,)^,3/4+5„/2/(2y^). (55) 

If the tensor force is attractive, the function / takes the 
form 

fA{x)^Cssmx + Cccosx, (56) 

with X = 2y/^ and where the specific linear combina- 
tion of sine and cosine factors is determined by the LO 
counterterm. In contrast, for a repulsive tensor force we 
have a regular and irregular solution. We have indeed 
that 

fnix) ^ Crcge^" + Ci,,e+^ . (57) 

where the solution is regular at the origin, while the 
e^^ solution diverges exponentially for r ^ 0. For the 
repulsive case we will always chose the regular solution 
at LO. 

In the attractive case, the convergence of the pertur- 
bative integrals is solely determined by the power-law 
behaviour of the reduced wave functions: the sine and 
cosine factors do not play any role on the finiteness of 
the perturbative corrections. Renormalizability requires 

5 1 

-n,>v+-, (58) 

which in turn implies that the first perturbative correc- 
tion to the C'2n.i counterterm scales as 

C2n.j ~ , (59) 

that is, the first correction to Cq enters at order Q~^^^, 
the C2 counterterm at order Q^, and so on. It is inter- 
esting to notice that (i) the power counting of the coun- 
terterms does not depend on the partial wave considered, 

(ii) the first perturbative correction to Co is ©(Q^^/^), 
which means that the infrared fixed point is stable, and 

(iii) the appearance of new counterterms is delayed in 
comparison with the non-singular potentials. 

In contrast, the analysis for repulsive tensor force is 
puzzling. The decreasing exponential behaviour of the 
wave function at short distances is able to render the 
perturbative integral finite, independently of the degree 
of singularity of he v-th order potential. This feature can 
be easily appreciated by considering the integral 

^^''^-/^^3T^-pM/?), (60) 

which is always finite. However, for a consistent power 
counting to emerge we need new counterterms at higher 
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orders. Two solutions are possible: (i) follow the power 
counting rules for attractive singular interaction, that is, 
C2n ~ Q(5"-l)/2, (ii) follow naive dimensional analysis, 
i.e. Weinberg counting, leading to C2n,i ^ Q^""*"^'. 

The first solution was proposed by Birse in Ref. [47| . 
It only considers as relevant to the power counting the 
power-law behaviour of the wave function, but not the 
sine, cosine or exponential factor. This expectation is 
consistent with the fact that the exponentials are just 
a sine or cosine function with an imaginary argument. 
However, it is also natural to expect that short distances 
do not play a significant role if the particles cannot ap- 
proach each other due to the repulsive potential barrier. 

This observation leads to the second solution, namely 
to follow naive dimensional analysis, which implicitly as- 
sumes that the repulsive channel can be treated as a per- 
turbation Q. For the particular case of the ^Pi channel, 
this is not an unreasonable prospect in view of the anal- 
ysis of perturbative OPE by Birse ^3], which suggests 
a critical momentum around 400 MeV above which the 
perturbative treatment will fail. In contrast, the criti- 
cal momentum for the attractive "^Pq channel is around 
200 MeV, requiring the non-perturbative treatment of 
the OPE tensor force if we want to describe the phase 
shifts up to /c ^ 2m7r. In addition, the critical mo- 
menta quoted above were obtained in the — limit. 
For finite pion masses we expect the critical momenta 
to rise, as the behaviour is strongly suppressed for 
TO^r > 1. In the repulsive case, the potential barrier pre- 
vents the two nucleons to explore in detail the m^r ~ 1 
region, where the singular nature of the intermediate 
range tensor force manifests, meaning that the critical 
momenta are expected to rise significantly with respect 
to the —J- limit. On the contrary, in the attractive 
case, the two nucleons are driven to the m^r ~ 1 region, 
and thus the critical momentum will only rise weakly. 

However, there is still the problem of what happens 
if the repulsive singular interaction is non-perturbative. 



Alternatively, we can try to analyze the scaling of the wave 
function at distances of the order of the breakdown scale of the 
theory, r ~ 0.5 fm. If the wave function behaves as r'+^ (instead 
of the singular r'^/* exp(— 2-^/03 /r)), naive dimensional expec- 
tations are fulfilled and the use of Weinberg counting is justi- 
fied. Unfortunately, for the particular case of the ^Pi channel 
the singular behaviour is achieved at distances of about 1 fm (as 
03 = 1.34 fm in this case), so the wave function scaling argument 
cannot be applied. The non-perturbative wave function scaling 
at the breakdown radius does not imply however the failure of 
the perturbative treatment of OPE. The reason is that the non- 
perturbative contribution (i.e. from r < 1 fm) to the phase shifts 
is particularly small, as can be explicitly checked from applying 
the methods of Rcf. 64] to this case, leading to 



drc 



,3 3/2 1/2 , 
' k rc 03 exp (- 



-)• 



The vanishing exponential ensures that, once the non- 
perturbative regime is achieved, the corresponding contributions 
to the total phase shift will be negligible. 



It is clear that the current understanding of the power 
counting of repulsive singular potentials is incomplete. 
Even though in the particular case of the ^Pi partial 
wave we can overcome the problem by invoking the per- 
turbative treatment of tensor OPE in this channel, the 
issue may eventually have important consequences for the 
power counting of coupled channels. For the moment we 
will simply assume that attractive and repulsive singular 
interactions follow the same power counting, unless we 
expect the singular potential to behave perturbatively. 
It may be possible that the extension of the methods of 
Long and Yang 57] from the '^Pq to the '^Pi channel will 
provide new insight into the power counting of repulsive 
singular interactions. 



C. Coupled Channels 

For analyzing the perturbative renormalization of chi- 
ral TPE in the coupled channel case, we employ the 
eigen representation of the phase shifts [7l|. In this 
parametrization the DWBA formulas take their simplest 
form. The eigen phase shifts are expanded according to 
the power counting as follows 

6^,{k; = £ Si]\k; r,) + 0{Q''--+') , (61) 



6p,ik; r,) = £ 4;.)(fc; r,) + ©(Q-^— +i) , (62) 



e,(fc;r,)= ^ e^;\k;r,) + 0{Q-'--+') , (63) 

and the LO phase shifts are computed non-perturbatively 
as explained in Appendix [X] The order v-th contribution 
to the eigen phase shifts is given by 



sin^4^) 



2^ 



A 



2/^ A-iy 



(fc)C(fc;^c) 



(64) 



{k)&hk-r,), (65) 



(-1)/ 



, 2^ A^^;\k)A-;\k) 
^ ^^^^^cot^-^^-cot^^) ^-^ 

(66) 

where the perturbative integrals laaj, Ipaj and /^g/jj are 
defined as 



dr 



(67) 
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with p, cr = a, /3; the subscripts a, b in the potential, 
which take the values a = j — 1 and h = j + 1, are used 
to denote the two angular momentum components of the 
potential (/ — j ± 1). In turn, we employ a, /3 to label 
the two asymptotic solutions at large distances, which 
behave as 



cos e 



(-1) 



(cot(5i/^ jj_i(fcr) - yj^ilkr)) , 
sin 4"^^ {cot 5^~l^ jj+i{kr) - yj+i{kr)) , 



(68) 



sme- -' {coi5^p.p jj-i{kr) - yj^i{kr)) , 



cose 



(-1) 



(cot 5^/^ j3+i{kr) - yj+i{kr)) , 



(69) 



where ,1^ . are normalization factors, which we define 

as fc^-i si;'' = -4i7^(fc) and fc^"+i = -^^^^'(fc). Fi- 

nally, for obtaining the perturbative corrections to the 
phase shifts in the nuclear bar representation 72] , which 
is the most widely used, we consider their relationship 
with the eigen phases 



Sij + 1 



2j 



sin(5ij - S2j) = 



tan 2€i 



tan 2ej ' 



sin 2ej 
sin2ei 



(70) 
(71) 

(72) 



to Up > 2, short distances are progressively more sup- 
pressed, which in turn implies that the paj integrals 
will eventually converge for large enough n. 

From the energy expansion of the reduced wave func- 
tions, Eqs. ([71)1 and ([75]) . we can see that the exact 
degree of divergence of the perturbative integral depends 
on the specific details of the channel and the phase shift 
under consideration. However, two simplifications are 
applicable for the particular case of the nuclear potential 
in chiral EFT. The first is that all the components of the 
subleading order potential diverge as l/r'^^", indepen- 
dently of the particular angular momentum subchannels 
which are involved. The second is that, for an iterated 
singular interaction at LO, the power law behaviour of 
the wave function at short distances does not depend on 
the angular momentum or the a or /3 nature of the scat- 
tering state: in fact we have s = Saa — Sba = — Sb/3 
and t = tact = ^ba — t&i3 = th/S- In this case, the situation 
is analogous to the uncoupled case and the behaviour of 
the perturbative integrals is given by 



2n,pa 



(re 



dr 



f,3-\-i' — 2s—nt 



(76) 



which diverges for 3 + — 2s — ?i< < 1. 

Independently of the previous simplifications, the es- 
sential point is that the perturbative integrals can be 
renormalized by the addition of a certain number of coun- 
terterms. As in the uncoupled channel case, the most 
explicit regularization method is to include nc(pa-) free 
parameters, that is 



E 

n=0 



+ 4:Ufc;^c), (77) 



and re-expand these relations according to the power 
counting. 

As in the uncoupled case, the renormalizability of the 
coupled channel phase shifts can be analyzed in terms of 
the fc^ expansion of the perturbative integrals 



n=0 



(rc)fc' 



(73) 



where only the first few terms in the expansion are di- 
vergent. The properties of the I'^^i paj terms are in turn 
related to the fc^ expansion of the reduced wave functions, 
which reads 



.Wk,pj{r)J 



t'o \'W2n.pj{r)) 



where U2n pj and W2n pj behave as 



U2n,p]{r) 

,W2„,pj(r) 



j,Sbp+ntbp 



(74) 



(75) 



at short enough distances. The specific values of sip and 
tip depend on the details of the LO potential. Owing 



where nc(pa) is the number of counterterms needed to 
regularize the pa integral. The total number of countert- 
erms is simply = nc(aa) + ^^c(a/3) + "c(/3/3)- 

The relation between the \^pa parameters and the 
counterterms can be obtained by postulating an explicit 
representation of the short range physics. A convenient 
representation is 

V^Sl(r-c) = ^'-''-^'^'^^^"^'^^ 



47rr? 



c(ll') 



E C(i,(r,)fc^"5(r-r,), (78) 

?i=0 

with fiir^) — ''"^^^P" and where nc{w) is the number 
of counterterms between the partial waves I and V . The 
relation between nc{ii') and nc(pa) will be discussed below. 
The momentum space representation of the the previous 
potential in the rc — > limit takes the form 



{p\v^t\p': 



E 

n=0 



U2n 



(79) 
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that is, they correspond with the usual representation 
of energy-dependent counterterms. If we have employed 
the short range normalization of Eq. (j32| , we obtain the 
relationship 

^2lapj = ^^C'2«!aa(^c)Wp.o(?'c)Ua.o(?'c) " 



^C'^;;^a6('^c)(up,o(?^c)Wa,o(?-c) + 



C^^nMi'^c)wp^o{rc)w^fl{rc) . (80) 

As can be seen, the previous representation of the short 
range physics requires n^^aa) = n^^aa), ^^(^ab) = n^^afi) 



and 



This is however not a universal 



feature: most regulators only require that the total num- 
ber of counterterms remains the same, and in general 
we should expect more counterterms in the lower partial 
waves than with the delta-shell regularization. 

A simplification that occurs in the coupled channel 
case is that the LO interaction is always singular, as a 
consequence of the tensor force being responsible of cou- 
pling channels with different angular momenta. More- 
over the singularity structure of the tensor force implies 
the existence of an attractive and a repulsive subchan- 
nel in all the coupled waves, that is, we do not need to 
discuss the two cases separately. The regular potential 
case will not be discussed either, but the results are the 
expected: if we consider the regular solutions, we end up 
with the standard Weinberg counting. On the contrary, 
in the particular case of the ^Si — ^Di partial wave, if we 
consider the irregular solution for the s-wave, the KSW 
counting appears. 



1. Singular Potential 

In coupled triplet channels, the tensor force piece of 
the OPE potential behaves at short distances (r — > 0) as 



2/.vL~l(r)-^±S, 



03 
J ^3 ' 



where V is a convenient matrix notation for the coupled 
channel potential, that is 



(82) 



with = a, b = j ± 1. The matrix elements of the 
tensor operator, S-', reads 



1 



- 2(j - 1) 6v<7C?~l) 

2J + 1 V6V7(J~1) -2(j + 2) ^ 



(83) 



Finally, is a length scale related to the strength of the 
tensor force. 



The interesting point for the regularization of the ten- 
sor force is that, at short distances, it can be diagonalized 
by means of the transformation 



2^Rj V. 



(-1) 




(84) 



which implies the existence of a repulsive and attractive 
eigenchannel. Due to the singularity, the tensor 

force overcomes the centrifugal barrier at short distances, 
and the reduced wave functions can be expressed as a sum 
of the attractive and repulsive solutions. In terms of the 
expansion of the reduced wave function, Eq. (|74t . the 
following behaviour is expected 



"2n,pj(r-) 



^ ^3/2+5n/2 



(85) 



where, depending on the sign of the tensor force, we can 
either have xa — 2-\/2a3/r, xr — l^f^a^Jr or xa — 
2-^/403/7", XR — l^pla^Jr. The functions gapj and gbpj 
follow the general pattern 

g{xA-,XR) = CssmxA + Cc cos + Crg''"^ , (86) 

where we are only taking the regular solution for the re- 
pulsive eigenchannel. The previous behaviour in turn im- 
plies that the repulsive solution plays no significant role 
in the renormalization of the perturbative phase shifts, 
as it only represents a negligible contribution to the LO 
wave function at short enough distances. 

From the previous form of the reduced wave functions, 
we have indeed the simplifications which we expected as 
a consequence of iterating a singular interaction at LO, 
namely that 



= Sb/3 

^b^g = 



= 3/4, 
5/2. 



(87) 
(88) 



If we compare with the s-wave regular solution of a reg- 
ular potential (s = 1, i = 2), the values above imply 
that the first perturbative correction to Cq appears half 
an order earlier than expected, while the C2n's happen 
(n — l)/2 orders later. The finiteness condition reads 



> V 



(89) 



which is independent of whether pa ~ aa, a/3 or /3/3. 
In terms of the scaling of the counterterms, the previous 
translates into 



C2n,ll' Q 



(5n-l)/2 



(90) 



that is, we obtain the expected scaling resulting from the 
iteration of a l/r^ potential. However, due to the exis- 
tence of the three different phase shifts, the number of 
counterterms triple with respect to the uncoupled chan- 
nel case. That is, at order Q'^l"^ we have a total of three 
counterterms, at order we reach six counterterms, and 
so on. 
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Even though the previous six counterterms guaran- 
tee the existence of the Tc — )■ limit of the NLO and 
N^LO scattering amphtudes, the existence of the repul- 
sive eigenchannel in tensor OPE means that finiteness 
may be probably assured with a smaller number of con- 
tact operators 0. The development of this argument 
requires to work in the basis for which the OPE tensor 
force is diagonal (i.e. the attractive-repulsive basis), in- 
stead of the usual partial wave basis. This is non-trivial 
from the technical side and the problem remains of which 
is the adequate power counting for the repulsive eigen- 
channel. However, if we simply accept the proposal that 
the repulsive components of the LO tensor force follow 
the same power counting as the attractive ones [i^l, the 
argument of finiteness becomes inconsequential from the 
power counting point of view and we end up with the 
aforementioned six counterterms at NLO/N^LO. This is 
the choice we are following for the coupled channels in 
this work. 

Alternatively, we can overcome the limitations in the 
understanding of the power counting of repulsive singular 
interactions if we suspect tensor OPE to be perturbative, 
as happened in the '^Pi partial wave, in which case we 
can use the Weinberg counting. There are two possibili- 
ties along this line: to work (i) in the attractive-repulsive 
basis or (ii) in the partial wave basis. The first case is 
interesting and leads to an exotic power counting 

E, but 

is also difficult to implement, as already commented in 
the previous paragraph, and will not be considered in the 
present work. The second case is much more straightfor- 
ward and physically compelling and will be analyzed in 
greater detail in Sect. IIVCI 

D. Convergence of the Perturbative Series 

The convergence of the perturbative series can be un- 
derstood in two different senses, namely with respect to 
the power counting expansion and with respect to the 
cut-off. The first kind of convergence is related with the 
EFT expansion parameter, while the second is linked to 
the interpretation of the cut-off within the EFT frame- 
work. 

We do not provide here definitive arguments, but 
rather educated guesses about these aspects of the the- 
ory. In particular, we can establish some bounds on the 



In particular, the two counterterms renormalizing the attractive 
eigenchannel may be enough. Of course, a formal proof will 
require to show that the divergences in the perturbative integrals 
are correlated. 

In particular we have C2n,AA ~ C2n,AR ~ 

Qi-5/4+2n jjj^j C2„,RR ~ Q2(i-l)+2n^ ^^leTe A and R refer to 
the attractive and repulsive eigenchannel respectively. It should 
be noted that the wave function in the repulsive eigenchannel is 
a mixture of the Uj, (' = i ~ 1) a^nd '^fc = i + 1) wave functions, 
and consequently it behaves as / = j — 1 in terms of the Weinberg 
counting. 



breakdown scale and the expansion parameter by consid- 
ering the related EFT expansion of the chiral potential 
in the first place. From the power series of the potential, 
see Eq. (jl9l) , it is apparent that the expansion parameter 
is formally a;o — Q/Aq (independently of which is the 
concrete value of Ao). Naively, we should also expect a 
power series in terms of xq for the scattering amplitude. 
However, we can argue on general grounds that the ex- 
pansion parameter of the amplitude will be bigger than 
that of the potential. That is, we have 

i/=—i ^ ^ ^ 

instead of Eq. (I20p . where the real expansion parameter 
is not xo = Q/Aq, but rather xi = Q/Ai, with Ai (< Aq) 
the true breakdown scale corresponding to the particular 
power counting under consideration. The departure from 
the naive expectation Ai = Ag is explained in terms of 
the iteration of the subleading pieces of the chiral poten- 
tial, which can spoil the convergence of the perturbative 
expansion of the scattering amplitude. 

The previous idea can be illustrated with the KSW 
counting [H, |4^ : in the singlet channel, the KSW break- 
down scale is determined by considering the effect of the 
iteration of OPE on the running of the counterterms [i^ , 
yielding Ai^s = Ann — 300 MeV. This value is much 
lower than the expected Aq ^ 0.5 GeV. The expansion is 
even less converging in the triplet channel: the analysis of 
perturbative OPE made by Birse [47j] indicates a break- 
down scale Ai^t — 100 MeV, in agreement with the results 
of Ref. E3- In this particular case we have a clear ex- 
ample in which the convergence is severely limited as a 
consequence of higher order perturbation theory . That 
is, the correct identification of the non-perturbative con- 
tributions to the interaction is an essential ingredient for 
a convergent EFT formulation. 

In the case of^interest for this work, the Nogga et 
al. counting f2l|, the breakdown scale can be esti- 
mated by the deconstruction method of Refs. (58l - [60| . 
in which the form of the short range interaction is de- 
termined by removing the pion exchange effects from 
the phenomenological phase shifts. This analysis sug- 
gests a breakdown scale of Ai „ ~ 270 MeV for the 
and ^Pi singlet channels [5^ [S^. In the triplet chan- 
nels, the deconstruction has only been performed in the 
uncoupled p- and d-wave triplets {^Pq, "^Pi and ^i?2), 
yielding Ai_t ~ 340 MeV [5a]. From the previous re- 
sults, we expect the approximate expansion parameters 
xi^s = 'tItt/Ai^s — 0.5 and xi,t = m-^j ~ 0.4 for the 
singlet and triplet channels respectively PI. 



In this regard, the KSW reformulation of Ref. [53 | claims to 
have solved the convergence problems of the triplet channel. 

It should be noticed that the deconstruction estimations for the 
breakdown scale are computed for a cut-off radius of rc = 0.1 fm. 
However, as is explained in the next paragraph, EFT calculations 
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The second type of convergence, that is, the conver- 
gence with respect to the cut-off, can also be stated in 
terms of the potential instead of the scattering amplitude. 
In this regard, the appearance of the 1 / r^'^'^ behaviour 
in the chiral expansion of the potential indicates that the 
potential is not convergent at short enough distances. 
If we consider the explicit expansion of the chiral po- 
tential at distances below the pion Compton wavelength 
{rriTrr < 1), that is 

1 °° 1 
FnnWoc-J:^, (92) 

this feature becomes apparent: the Taylor expansion does 
not converge if r < i?o ~ l/^o- Within the framework 
of a non-perturbative treatment of the chiral potential, 
this argument requires the coordinate space cut-off to be 
bigger than the breakdown radius, Vc > Ro, which we 
naively expect to lie in the Rq ~ 0.1 — 0.5 fm region. 
Otherwise the scattering amplitudes will start to diverge 
at high enough orders, independently of the number of 
counterterms included in the computations. This in turn 
sets a limit on the value of the cut-off in perturbative 
calculations. Again, as a consequence of subleading order 
iterations, the bound of the cut-off will become softer, 

R\ > i?o- 

From this observation, a natural interpretation of the 
cut-off arises within the context of approximation theory, 
as originally advocated by Stevenson jTSj] . While the full 
EFT scattering amplitude (computed at infinite order) is 
independent of the cut-off in the regions where the poten- 
tial converges as a consequence of containing an infinite 
number of counterterms, the truncated EFT scattering 
amplitude (computed at finite order) is cut-off depen- 
dent. However, as far as we have included the necessary 
counterterms guaranteeing renormalizability 0, we are 
free to chose whatever value of the cut-off that provides 
the better convergence properties towards the full result 
and realizes the particular power counting under consid- 
eration. In this sense, the cut-off is just a parameter con- 
trolling the con verg ence of the theory, as recently stated 
by Beane et al. |54| . 



will not converge below a certain value of the coordinate space 
cut-off probably around 0.5 fm. Although results computed be- 
low that point are not to be trusted, the truth is that the de- 
construction estimations are stable up to Tc = 0.8 fm, which is 
within the values of the cut-off which we regard as acceptable 
in this work. In this regard, we expect the final value of the 
expansion parameter to be similar to the estimations based on 

Refs. [H-iai. 

This condition is required for the difference between the full re- 
sult and the calculation at i^-th order to formally scale as Q"^^^, 
see also the related comments at the end of Sect. Ill CI 



E. Power Counting 

The power counting arising from the perturbative 
treatment of chiral TPE is summarized in TableHl The to- 
tal number of counterterms in the Nogga et al. counting 
at NLO and N^LO {uc = 21) is certainly bigger than the 
corresponding one in the dimensional Weinberg count- 
ing (uc = 9). The counterterm pattern is fundamentally 
very similar to the power counting obtained by Birse us- 
ing RGA '47']. However, there are two minor differences 
with respect to the results of Ref. 01 that are worth 
commenting: (i) the counting of the ^5*1 partial wave, 
and (ii) the size of the C2n operators. 

In the RGA of Ref. [l^l, three counterterms (instead 
of two) are predicted for the '^Si channel. As explained 
in Ref. [ssj . this discrepancy arises as a consequence of 
the naive extension of the idea of trivial and non-trivial 
fixed points of the RG equations [g^, 113] , which is only 
valid for regular long-range potentials, to the singular 
interaction case in Ref. |47| . However, from the point 
of view of perturbative renormalizability, the different 
fixed points stem from the power-law behaviour of the LO 
wave functions at short distances. In this regard, regular 
long-range potentials require the existence of two fixed 
points as a consequence of the existence of a regular and 
irregular solution. That is, if we have a regular solution 
the matrix element of the counterterm is 

(^c)rcg = C(r,) J dr^S{r r,) ^ C(r,) rf , 

(93) 

which requires C{rc) to scale as r~^' (i.e. Q^') for having 
a non-trivial effect. On the other hand, for the irregular 
solution we have instead 

(^c)irr^C(r,)r-2'-2, (94) 

explaining the Q^^'"^ scaling. However, in the partic- 
ular case of attractive singular potentials, all the short 
distance solutions of the wave function show the same 
power law behaviour, and therefore there is only one 
type of power counting in this case, which is given by 
the Q"^/^ scaling as can be trivially checked by the pre- 
vious procedure. 

The second difference lies in the order of the energy- 
dependent operators, C2n, which is higher in the present 
formulation (C2„ - Q(5"-i)/2) than in the RGA of 
Ref. [13 (C2„ - Q(4n-i)/2)^ A possible explanation 
is that the expansion of the reduced wave function 
for the attractive triplets generates a stronger power law 
short range suppression than the corresponding one for a 
regular potential, r^/^ for each power of fc^ instead of r^, 
shifting the C2n operators to higher orders by a factor of 
Q"/^. Taking into account the previous effect, which may 
have been overlooked in the analysis of Ref. [43| as a re- 
sult of the presence of two cut-offs, the C2n operators are 
demoted from Q(4n-i)/2 Q(5n-i)/2^ reconciling RGA 
and perturbative renormalizability. However, the recon- 
ciliation is also possible if we consider energy-dependent 
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contributions to the chiral potential, that is 

h,2m 



(95) 



where 2m, < k — \/ /i -Eiab is the center-of-mass mo- 
mentum and -Blab the laboratory energy. This contribu- 
tions are not forbidden by power counting, but do not 
occur in the present formulation which employ energy- 
independent potentialsEl- If they are taken into account, 
the counterterms should be promoted from (2(5n-i)/2 -(-q 
Q(4n-i)/2 perturbative framework, recovering the 

initial result from Ref. Finally, a third explanation 
is provided by the possibility that the r^/^ suppression 
related to each subtraction may be contaminated by the 
low energy scale Aqpe- In such a case, the missing scale 
would imply the reinterpretation of r^l"^ as \/r Aqpe, 
thus enhancing the C2n counterterm by Q^"/^. 

In principle, the discussion is merely academical as it 
only affects the counting of counterterms beyond the or- 
der considered in the present work. However, at N'^LO 
(Q*) taking one option or the other can substantially 
change the counting: if we assume the scaling C2n ~ 
q(5"-i)/2^ the third counterterm in triplet channels does 
not appear until order Q^/^, between N^LO and N'^LO. 



On the contrary, if we follow Ref. [47[ , the triplet chan- 
nel C4 counterterms will enter at Q*, increasing the total 
number of free parameters from 27 to 35 at N'^LO. As 
the energy-independent representation of the finite range 
potentials is the preferred one (and as far as the Aqpe 
contamination hypothesis has not been checked), we ad- 
vocate for the first option. In this regard, the recent 
power counting analysis of the ^Pq channel by Long and 
Yang 57] also prefers the scaling (instead of Q'^^^) 
for the C2 operator of the attractive triplet, in agreement 
with our observations. 



IV. THE P- AND D-WAVE PHASE SHIFTS 

In this section we compute the nucleon-nucleon p- and 
d-wave phase shifts within the perturbative framework 
developed in the previous section, that is, we treat the 
LO potential non-perturbatively and include the sublead- 
ing order corrections as perturbations. For the countert- 
erms, we follow the results of Table HI The details of the 
calculations are described in the following paragraphs. 

The exact form of the chiral nucleon-nucleon poten- 
tial in coordinate space is taken from Ref. [T^I- We 
use the parameters gA = 1.26, f-^ = 92.4 MeV, = 
138.03MeV, and dig = -0.97 GeV^ The effect of the 
Goldberger-Treiman discrepancy (dis) in the chiral po- 
tential can be effectively taken into account by chang- 
ing the value of gA from 1.26 to 1.29 in the LO piece 



With the exception of the relativistic corrections to one pion 
exchange, for which the energy-dependent representation is cho- 



of the potential, see Refs. [l^ [3| for details. For the 
chiral couplings we use the values ci = — 0.81 GeV""'^, 
C3 = —3.40 GeV~^ and C4 = 3.40 GeV~^, which are com- 
patible with the determination of Ref. 75] . We include 
the relativistic corrections to OPE at NLO and the re- 
coil corrections to TPE at N^LO. Following Ref. [lil, the 
relativistic corrections to OPE are included in an energy- 
dependent manner, that is 



(2) 



OPE 



(r)=- 



2M2, 'OPEV 



(96) 



where k is the center of mass momentum and A/^r = 
938 MeV is the nucleon mass. It should be noted that the 
s-wave results of Ref. [HI did not include the relativistic 
correction to OPE. However, explicit calculations show 
that the effects of the relativistic corrections are fairly 
small: if we include the contribution from Eq. ([M)) . the 
s-wave phase shifts do not change appreciably. The same 
observation applies if we remove the recoil 1/Mjv TPE 
corrections from the results of Ref. [11]. In fact, if we 
consider the deconstruction estimation for the breakdown 
scale of the theory, that is, Ai^^ ^ 300 MeV and Ai_t ~ 
350 MeV, we have the approximate relationship 



Q_ 

Mn 



(97) 



supporting the original convention of Weinberg for the 
Q/Mn terms 

The LO phase shifts are computed non-perturbatively 
by iterating the OPE potential to all orders and by 
adding the counterterms required to achieve cut-off in- 
dependence at this order [21|. This condition requires 
the inclusion of a counterterm in the ^Sq, '^S'l — ^Di, 
"^Pq, ^P2 — and ^D2 partial waves at LO. In this 
work the counterterms are included in an implicit man- 
ner by fixing the value of the scattering length in the LO 
phase shifts for the previous channels. The technical de- 
tails are explained in Appendix lAl Following Ref. fEE], we 
employ the cut-off window Vc ~ 0.6 — 0.9 fm for regular- 
izing the phase shifts. In addition, we show the results 
for Tc — 0.3 fm for comparison, which also serve as an 
informal check of the hard cut-off limit of the scattering 
amplitudes. However, contrary to what happens in the 
s- waves |55j] , the variation of the phase shifts in the previ- 
ous cut-off window can be hardly interpreted as the error 
band of the results at a given order. The reason is that 
the cut-off dependence of the p- and d-wave phase shifts 
is fairly small, sometimes negligible, for cut-off radii be- 
low rc — 1.2 fm. 

The subleading order corrections to the phase shifts 
are computed by making use of the perturbative for- 
malism developed in the previous section. This formal- 
ism requires a total of 21 counterterms at NLO/N^LO, 
distributed among the partial waves according to the 
results of Table HI These counterterms, or equivalently 
the free parameters which modify the perturbative in- 
tegrals, see Eqs. (|1T|) and ([771) . are determined by fit- 
ting the NLO and N^LO perturbative phase shifts in 
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FIG. 3. (Color online) Phase shifts (nuclear bar) for the ^Pi, ^Po, and ^P^ — F2 channels with the LO piece of the 
chiral potential (OPE) fully iterated and the NLO and N^LO pieces (chiral TPE) treated perturbatively. A LO counterterm in 
included in the attractive triplets (^Po, '''P2 — ^F2) to remove the uncontrolled cut-off dependence. The contact operators are 
used to fix the value of the scattering lengths to 0,3 — — 2.7fm^ and asp^ = — 0.04fm^. The bands are generated by varying 
the cut-off radius in the region = 0.6 — 0.9 fm and they are formally a higher order effect. However, as a result of the mild 
cut-off dependence of the results, the full uncertainty of the calculation at a given order is expected to be much higher than the 
cut-off variation bands. The light blue band corresponds to the N^LO results of Ref. [l^ in the standard Weinberg counting. 
The dashed dark blue line represents the N'^LO results for Vc ~ 0.3 fm. 



the k = 100 — 200 MeV region to the Nijmegen II phase 
shifts [76], which are in turn equivalent to the Nijmegen 
PWA [77]. 

It should be noticed that the perturbative regulariza- 
tion techniques employed in this work are specifically cho- 
sen to identify divergences in the amplitudes rather than 
to optimize the phenomenology. In fact, the regulator 
we employ can be considered to be the coordinate space 
equivalent of the sharp regulator in momentum space, 
which is known to be suboptimal from the phenomeno- 
logical point of view. Even with this proviso, the phase 
shifts we obtain are better than the corresponding ones 



in the Weinberg scheme at the same order. However, the 
amplitudes are clearly amenable to improvements. 



A. P- waves 

The results for the p-wave phase shifts are shown in 
Fig. ([3]). The renormalization of the LO phase shifts 
requires the inclusion of a counterterm in each of the 
attractive triplets. Therefore we fix the value of the 
scattering length in the ^Pq and ^P2 — '^^2 channels to 
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FIG. 4. (Color online) Phase shifts (nuclear bar) for the ^D2, '^D2 and ^Ds — ^Gs channels, with the LO piece of the chiral 
potential (OPE) fully iterated and the NLO and N'^LO pieces (chiral TPE) treated perturbatively. A LO counterterm is 
included in the ^D2 channel, as suggested in Ref. [2l|, which is used to fix the scattering length to the value asjj^ = — 7.4fm^. 
The bands are generated by varying the cut-off radius in the range = 0.6 — 0.9 fm (see the related comments about their 
interpretation in the previous figure and in the main text). The light blue band corresponds to the N^LO results of Ref. [l^ 
in the standard Weinberg counting. The dashed dark blue line represents the N^LO results for rc = 0.3 fm. 



asp^ = -2.7fm^ and aap^ = -0.04fm^. The previ- 
ous values differ from the scattering lengths that can 
be obtained from the Nijmegen II potential, namely 
aspjj Nijjjj = — 2.468 fm'^ and asp^ Ny^n = — 0.2844 fm'^ ac- 
cording to Ref. The LO deviations of the "^Pq and 
— scattering lengths from the Nijmegen II values 
are however inconsequential: the LO values should only 
be accurate up to higher order corrections and there is 
no need to use the exact value in the LO calculation. 
The important point for the LO scattering lengths is to 
provide a good starting point for the power counting ex- 
pansion. In addition, the mismatch between the periph- 
eral effective range parameters of the Nijmegen poten- 



tials and the preferred values of the chiral potential is 
well known [13, H, iH, [H . 

In general, the description of the p- waves is better 
than in the Weinberg counting, although this is a natural 
prospect once we take into account the extra countert- 
erms. In this regard, the ^Pi and ^Pi wave are interest- 
ing in the sense that they only contain the counterterms 
prescribed by Weinberg's original counting. Therefore 
they provide a more direct comparison with the Wein- 
berg scheme, and show that there is no significant draw- 
back to treating TPE perturbatively in these waves. In 
fact, the results are slightly better than in the standard 
Weinberg counting. 
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The phase shifts are nicely reproduced up to around a 
center-of-mass momentum of fccm = 300 MeV, ahhough 
there are signs of convergence up to kcm = 350 MeV, a fig- 
ure compatible with deconstruction '58"] . Contrary to the 
s-wave case, the value of the p-wave phase shifts in the 
chosen cut-off window does not differ much from those 
obtained when the cut-off is removed. Consequently 
we expect the deconstruction estimations, which are ob- 
tained for Tc = 0.1 fm to work better in the p- waves than 
in the s- waves. 



B. D-waves 



Partial wave 


LO 


NLO 


N^LO 


N^LO 




1 


2 


3 


3 




1 


2 


2 


3 


All 


6 


20 


21 


26 



TABLE II. Partial (and total) number of counterterms in 
the ^P2 - ^F2 and ^Ds - ^Gs partial waves at LO (Q~^), 
NLO (Q^), N^LO (g^) and N^LO (g*) assuming that (i) 
the iteration of the OPE potential is restricted to the lower 
partial wave of the coupled channel ("^^2, '^D^) and (ii) energy- 
dependent counterterms. As in the previous case, for the 
N^LO counterterms see the related discussion of Sect llll El 



In Fig. Q we show the results for the d-wave phase 
shifts. In principle, the complete renormalization of the 
LO d-wave scattering amplitudes requires the inclusion 
of a counterterm in the two attractive triplets: the 
and the ^D^ — "^G^ partial waves. However, we follow 
here the original suggestion of Ref. [2l|, in which the 
counterterm is only (optionally) incorporated in the 
channel. The reason is that regulator dependence in the 
•^Ds — ^Gs channel is not apparent until relatively high 
values of the cut-off [2l|. Thus, we can safely obviate the 
counterterm, at least in the LO calculation. On the other 
hand, the full iteration of OPE in the ^D^ — '^Gs channel 
will require the inclusion of six additional counterterms, 
largely reducing the predictive power of the theory. 

In the '^D2 channel we fix the value of the scattering 
length to asjj^ = — 7.4fm^, which basically coincides with 
the Nijmegen II value as = -7.405 fm^ [78]. In the 
^D^ — ^Ga channel, we impose regular boundary condi- 
tions for the LO amplitudes. This procedure is of course 
non-renormalizable, as commented in the previous para- 
graph, although the LO cut-off dependence does not ap- 
pear until we reach Vc — 0.2 fm. However, the NLO and 
N^LO cut-off dependence becomes manifest at higher val- 
ues of the cut-off, around Vc = 0.5 fm in the N^LO case. 
That is, the present treatment of the ^D^ — ^G^ phase 
shifts is inconsistent. This is also the reason why we do 
not show the Tc = 0.3 fm results for the ^D^ — ^Ga partial 
wave in Fig. (jH). 

If we do not include the LO counterterm, the correct 
renormalization of the '^D^ — '^Ga channel should require 
the OPE potential to be of order instead of Q^^. 
Consequently, the consistent iteration of OPE requires 
the evaluation of (v + l)-th order perturbation theory to 
calculate the Q'^ scattering amplitude, a procedure which 
would generate renormalizable results. On the contrary, 
the inconsistent full iteration of OPE at LO causes the 
amplitudes to eventually diverge as 



dr ■ 



r-3/2 

,3+1^ 



1 



1/2+1/ 



(98) 



once the LO wave function attains the expected short 
range behaviour for a solution of a potential. Of 
course, the implicit assumption behind the present cal- 
culation of the — '^Gs phase shifts is that the full- 



iterated (or inconsistent) results are not going to be sub- 
stantially different from the partially-iterated (or consis- 
tent) ones, provided the cut-off range is soft enough. For 
the range of cut-offs we employ {vc = 0.6 — 0.9 fm) the 
divergent regime has still not been reached, meaning that 
the current procedure probably represents a fairly good 
approximation of a consistent result. In any case, the 
elucidation of this aspect definitively calls for a pertur- 
bative reanalysis of the '^£'3 — ^G^ partial wave in order 
to check this assumption. 

It should be commented that the d-wave (and to a 
lesser extend the f-wave) scattering amplitudes have been 
traditionally a problem in nuclear EFT, especially at 
N^LO: the excessive attractiveness of the central con- 
tribution to the chiral potential at this order generates 
phase shifts which do not reproduce the results of the par- 
tial wave analyses, as first noticed by Kaiser et al. [7l,^3 ■ 
In particular this poses an issue with respect to the con- 
vergence of nuclear EFT in these waves: if the N^LO re- 
sults are still considerably different from the PWA ones, 
then we should still expect large corrections from higher 
order contributions and, consequently, a slow conver- 
gence rate, rendering the application of nuclear EFT im- 
practical. In this regard, Entem and Machleidt [l^ no- 
ticed that the inclusion of a counterterm in each of the 
d-waves in the N^LO potential effectively solves the prob- 
lem, as the short range operators provide the necessary 
repulsion for compensating the excess of attraction of the 
long range potential. Curiously, this is not very dissimilar 
to the situation we encounter within the present power 
counting scheme. A different solution which is more com- 
patible with the application of the Weinberg counting is 
the use of spectral function regularization, as advocated 
by Epelbaum et al. [3, [l^. In spectral function regu- 
larization a second cut-off is introduced to regularize the 
momentum in the pion loops, a procedure which greatly 
reduces the strength of the N^LO potential, effectively 
solving the convergence problem without the addition of 
new counterterms. Of course, the most obvious solution 
is the use of softer cut-offs, as happens for example with 
the Nijmegen xPWA's [zllIII- 

As we can see, the previous problem manifests mod- 
erately in the d-wave phase shifts of Fig. (U). In partic- 
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ular, it only affects the ^i?2 and phases, although 
in a milder form than in Ref. [795. any case, the 
excess of attraction in the and ^D^ channels for 
the — 0.6 — 0.9 fm cut-off window is still under con- 
trol, and we do not expect to have convergence problems 
for fccm < 300 MeV at N^LO, when the first countcrt- 
erm is added. However, it is clear that a softer cut-off 
window will probably improve the convergence rate of 
these two partial waves without worsening the other p- 
and d-wave phase shifts. Finally, the "^1)2 partial wave 
is well-described and shows signs of convergence up to 
^cm = 350 MeV, as expected from deconstruction [58] . 



C. Improved Power Counting Scheme for the 
Coupled Channels 

A particularly problematic feature of iterating OPE in 
the coupled triplet channels is the associated requirement 
of including a total of six counterterms. For the ^Si — ^Di 
channel, this condition is not a significant drawback. 
However, the counterterm proliferation caused by the 
full iteration of OPE starts to become worrisome in the 
— channel and excessive in the ^D^ — ^G^ partial 
wave. The reason for the six counterterm figure is that 
the iteration of OPE in the coupled channels assumes the 
same short-range behaviour for the two- linearly indepen- 
dent LO wave functions which correspond to the a and 
(3 scattering states. This justifies the necessity of two 
counterterms per phase: the scaling of counterterms is 
identical independently of whether they only affect the 
i^Dj) partial wave, the '^i^2 i^Gs), or the admixture 
between them. However, it is natural to expect that the 
r^/^ wave function behaviour will be first achieved in the 
a (rather than in the /3) scattering state. 

Along the previous lines we analyze here an alternative 
power counting scheme for the peripheral coupled triplets 
which (i) reduces the number of counterterms without 
significantly worsening the phenomenological description 
of the phase shifts in the •^P2 — ^F2 channel and (ii) limits 
the counterterm proliferation if we decide to iterate OPE 
in the '^D^ — ^G^ partial wave. The idea is to only iter- 
ate tensor OPE in the a scattering state of the coupled 
channel, that is, in the lower partial wave. For the partic- 
ular case of the ^P2 — ^F2 channel what we mean is that, 
while we have naively expected the full OPE potential to 
be promoted to order Q^^, that is 



' fp 



it is much better to do the following promotion 



y 'D'D pf 



(-1) y(0) 



(0) 



ff 



pp 
fp 



Pi 



V, 



(0) 



(99) 



(100) 



in which only the diagonal p-wave piece of the poten- 
tial is iterated. In this scheme, we have one and two 



(p- to p-wave) counterterms at orders and re- 
spectively, while the first counterterm mixing the p- and 
f- waves enters at order Q^+'^Z'*. Taking into account the 
smallness of the £2 mixing parameter and the S3 phase 
shift, the previous modification does not seem an unrea- 
sonable prospect. An advantage of this scheme is that 
the size of the matrix elements of the two-pion exchange 
potential are reduced in the pf and // channels, thus 
decreasing the size of the subleading contributions with 
respect to the naive implementation of the counting of 
Ref. [2l|. For the ■^£'3 — ^G^ partial wave we propose the 
following promotion 




(-1) T/(0) 



fo) 



(101) 



which implies one counterterm at LO for the dd subchan- 
nel, and two at NLO and N^LO. The first counterterm 
mixing the d- and g-wave subchannels appears at order 
q3+3/4^ that is, between N^LO and N^LO. The power 
counting resulting from these modifications is summa- 
rized in Table HIl 

The power counting proposed here can be explained 
on the basis of a variation of a well-known argument [2l|, 
'57'! that explains the perturbative character of OPE in 
peripheral waves. We analyze the behaviour of the partial 
wave projection of the OPE potential at low momenta, 
which is given by 



(p|VbpE,ii'|p'> 



P P 



(102) 



where Aqpe is the scale governing the strength of OPE 
potential, I and I' are the angular momenta of the partial 
waves and cw a geometrical factor that approximately 
behaves as 1/(Z!/'!) |21|]. If Aqpe is considered to be a 
low energy scale, the order of the OPE potential will be 
(instead of the naive Q"). However, owing to the far 
from perfect separation of scales, the geometrical factor 
Cw is eventually able to mimic the effect of positive pow- 
ers of Q/Aq (or more properly Q/Ai, see Sect IIII D] for 
details) Ej. In the nuclear EFT proposed in the present 
work, we expect a expansion parameter of the order of 
0.4 for the triplet, which means that for I oi I' > l„it 
(with Zcrit = 2 or 3) the angular momentum suppression 
is equivalent to the demotion of OPE by at least one 
order. For higher partial waves the actual demotion of 
OPE will be even higher, proportional to several orders 
in the chiral expansion, thus explaining the fact that first 
order perturbation theory for OPE works better the more 
peripheral the partial wave [79| . For keeping matters sim- 
ple, we have only considered the demotion of OPE of by 
at most one order for ^ > 3. 



There exists the possibihty that a similar geometrical suppres- 
sion is taking place for repulsive tensor forces when treated non- 
perturbatively, reconciUng power counting with the naive expec- 
tation that short range interactions should be less important in 
the repulsive case compared to the attractive one. 
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FIG. 5. (Color online) Phase shifts (nuclear bar) for the P2 — ^F2 and "^Dg — "^Gs coupled channels in a power counting scheme 
in which the OPE potential is fully iterated only in the lower partial wave of the coupled channel, while the other contributions 
are treated as perturbations. The LO counterterm are employed to fix the value of the scattering length to aap^ — — 0.04fm'^ 
and a3£)3 = — 0.085 fm^. The light blue band corresponds to the N'^LO results of Ref. [31 in the standard Weinberg counting. 
The dashed dark blue line represents the N^LO results for rc — 0.3 fm. 



Of course, the consistent implementation of the pre- 
vious power counting scheme requires the evaluation of 
perturbation theory up to third (fourth) order at NLO 
(N^LO). However, owing to the computational difficul- 
ties related to the evaluation of the perturbative series 
at high orders, we will fully iterate the OPE poten- 
tial as in previous cases. Although this is inconsistent 
from the power counting point of view, and probably 
even non-renormalizable, we expect this simplified cal- 
culation to be a good approximation of the true results 
in this improved counting scheme. The regularization of 
the subleading order amplitudes is performed by assum- 
ing a delta-shell parametrization of the counterterms, see 
Eq. ([75|. The corresponding A-parameters to be added 



in the perturbative integrals can be obtained by means 
of Eq. ([50]) . The counterterms are then fitted to the 
Nijmegen II phase shifts in the fccm = 100 - 200 MeV 
region. The amplitudes generated by this particular 
regularization procedure seem to be fairly cut-off in- 
dependent: no sign of divergent behaviour has been 
found for Vc > 0.3 fm. Numerical limitations have pre- 
vented the exploration of the harder cut-off region below 
0.3 fm. Consequently we are unable to determine the 
eventual renormalizability of the regularization we have 
employed F^ . 



Conversely, the direct inclusion of A-parameters only in the I3 
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The results are shown in Fig. ([5]). The LO phase shifts 
have been regularized by fixing the scattering length 
to aap^ = -0.04 fm^ and a-^/jg = -0.085 fm^. In the 
— ■^i^2 channel the improved power counting scheme 
yields similar results to the original one in Fig. ([3]), but 
employing half the number of parameters. The results 
for the '^I?3 — ^Ga have also improved, as expected from 
the inclusion of two new counterterms. The convergence 
properties are similar to other attractive triplets, and the 
Tc ~ 0.3 fm results do not significantly differ from the 
r„ = 0.6 — 0.9 fm ones. 



V. DISCUSSION AND CONCLUSIONS 

In the present work we have explored in detail the ex- 
tension of the power counting proposal of Nogga, Tim- 
mermans and van Kolck "21] to the subleading orders of 
the chiral expansion up to N'^LO (Q^)- We have chosen 
a perturbative regularization framework which facilitates 
the identification of divergent contributions and therefore 
the subsequent identification of the correct power count- 
ing for the contact operators. In particular, the present 
scheme fulfills the demanding constraints of renormaliz- 
ability and power counting without neglecting the phe- 
nomenological aspects: the results are in fact as good as, 
if not better than the corresponding ones in the Weinberg 
counting at the same order. 

We have determined the power counting of the con- 
tact operators, resulting in a total of 21 counterterms for 
the S-, p- and d-waves at NLO and N^LO, that is, twelve 
more than in Weinberg dimensional counting at the same 
order. The distribution of the contact operators largely 
agrees with the related RGA by Birse |43|, but we also 
include several clarifications and improvements. In the 
^Pi channel we explain the power counting disagreement 
between RGA ^3] and the original exploration by Nogga 
et al. [2l[. We also propose modifications which enhance 
the consistency of the approach and improve the coun- 
terterm distribution in the '^P2 — ^F2 and '^Ds — ^G^ chan- 
nels. However, a serious assessment of power counting in 
the previous cases requires a perturbative reanalysis of 
the OPE potential up to fourth order in the perturbative 
series. 

In principle, the perturbative treatment of the sublead- 
ing pieces of the chiral interaction raises the problem of 
a serious departure with respect to the traditional ap- 
proach to nuclear physics, namely the idea of using a po- 
tential as the elementary computational building block in 
few and many body calculations. However, the divorce 
is merely apparent. Perturbative renormalization should 
be understood as an analysis tool of the formal aspects 
of the theory, rather than as the adequate framework for 



final calculations. The point is not to explicitly treat the 
subleading orders of the chiral two-body forces pertur- 
batively, but the fact that they only represent a small 
contribution to physical observables. 

In this respect, a sensible approach is to construct chi- 
ral two- and three-nucleon potentials that can be shown 
(i) to contain the subleading pieces of the interaction ap- 
proximately as perturbations, therefore complying with 
power counting expectations, and (ii) to be cut-ofl[ inde- 
pendent modulo higher order uncertainties in a reason- 
able range of cut-offs. Non-perturbative methods may 
fulfill these conditions provided an adequate regulator 
function and cut-off window are employed, as proposed 
in Ref. [6l|. In particular, the perturbative character 
of chiral TPE can always be checked a posteriori. If 
we consider t he p ioneering work of Ray, Ordoiiez and 
van Kolck [111, Il2j | as the first generation of chiral poten- 
tials, and the phenomenologically successful potentials 
of Entem and Machleidt il7j and Epelbaum, Glockle and 
Meifiner [g^] as the second generation, what we need is a 
third generation of chiral potentials which overcomes the 
theoretical deficiencies of previous attempts. This task 
will surely require much less work than the remaking of 
all of nuclear physics in a perturbative fashion. 

The prospects for the feasibility of this approach are 
indeed good. In this regard, there exists the exciting pos- 
sibility of reinterpreting the Nijmegen xPWA99 |2J] as an 
EFT calculation in the present power counting scheme. 
In the chiral PWA of Ref. [7J| , the authors regularize the 



chiral N'^LO potential with a total of 23 (22) boundary 
condition at the cut-off radius Vc = 1-4 (1.8) fra, achieving 
a x^/d-o.f. ~ 1 in the -Biab < 350 MeV region[3- Taking 
into account the size of the cut-off radius they use, chiral 
TPE is almost guaranteed to be a perturbation, espe- 
cially in view of the perturbative analysis of Ref. [82l | . In 
the power counting advocated in this work, we employ a 
total of 21 parameters for the s-, p- and d-waves, a figure 
suspiciously similar to Ref. [74]. However, explicit calcu- 
lations would be required to check whether the xPWA99 
can be considered a realization of the counting proposal 
of Nogga et al. [21i] . 
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(Isjy^) and eventually the 1^2 perturbative integrals generate 
divergences in all the integrals that do not contain at least two 
A-parameters. 



^® Note that we are not taking into account the more recent 
XPWAOS of Ref. [8l|. which claims to describe the nucleon- 
nucleon scattering data up to Eiab = 500 MeV with a total of 
33 boundary conditions at Tc = 1.6 fm. The reason is that en- 
ergy range of Ref. 81] extends well beyond the expected range of 
applicability of nuclear EFT. 
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Appendix A: Renormalization of the Leading Order 
Phase Shifts 



In this appendix we explain the non-perturbative 
renormahzation of the LO scattering amphtudes and 
wave functions. The appendix can be considered a re- 
view of the coordinate space renormahzation techniques 
of Refs. [H-^. We consider first the uncoupled channel 
case and then move to the coupled channel case. 



1. Uncoupled Channels 



We start by defining the reduced Schrodinger equation 



for the LO reduced wave function 



2fi 



1(1 + 1) 



'■k,l 



(Al) 



which we consider to be valid from from Tc to cx), where Tc 
is the cut-off radius. With k and I we refer to the center 
of mass momentum and the orbital angular momentum 
respectively, and T^(~^)(r) is the LO potential in config- 
uration space. The asymptotic (r — > oo) normalization 
of the Uk wave function is given by 



where 5\^^\k;rc) is the LO phase shifts. 

The regularization and calculation of the phase shifts 
depend on whether we include a counterterm at LO or 
not. We will only consider in detail the first case, as the 
second can be handled by standard means. The addition 
of a counterterm is equivalent to the condition of fixing 
the scattering length in the zero energy solution. For this, 
we set the asymptotic form of the zero energy solution to 



(r) ^ fc' (cot 5[-^^ ji [kr) ~ yi (kr) ) , (A2) 



where a| is the LO value of the scattering length, 
which we fix to the desired value. We integrate the zero 
energy solution downwards, from oo to Tc in Eq. (jAip . 
Then, we assume the logarithmic derivatives of the zero 
and finite energy solutions to be equal at r = rc 
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which effectively provides the initial integration condi- 
tions for obtaining the finite energy solution at arbitrary 
radius r. Finally we integrate upwards and obtain the 



phase shifts by comparing the wave function at large radii 
with the asymptotic form of Eq. (|A2I) . 



However, for the previous wave functions to be useful 
in determining the divergence of the perturbative inte- 
grals we need to change the normalization. In the nor- 
malization defined by Eq. (|A2p . the reduced wave func- 
tion u^, is not necessarily energy independent at the 
cut-off radius, making more difficult the identification of 
the required counterterms. We therefore define a new 
normalization of the wave function which is energy inde- 
pendent &t r — Tc 



Uk,i{r) ^ A\ ^\k;rc)ukj{r) . 



(A5) 



The normalization factor a\ {k; Tc) can be uniquely 
determined by requiring the additional condition 
A[~^\0;rc) = 1, in which case we have 
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2. Coupled Channels 



For the coupled channel case the reduced Schrodinger 
equation reads 



2fi 



+2^,v},'>{r)w)^ 



(A7) 



-w 



2/x 



+2^.vL-»(r)u'-»=k^4-•\ 

(A8) 



where j is the total angular momentum and u (w) the 
I — j ~1 {I = j + 1) wave function. For the LO potential 
we employ the a (b) subscripts to denote the I = j — 1 
(l = J -l- 1) subchannel in the potential matrix. Of course, 
we assume the Schrodinger equation to be valid only from 
r = Tc to oo. 

For the coupled channel case we can define two linearly 
independent solutions for r — > oo, which are the a and /3 
scattering states. In the Blatt-Biedenharn parametriza- 
tion of the phase shifts 71] their asymptotic behaviour 
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is given by 



X (cot5i/)j,_i(fcr)-y,_i(fcr)), (A9) 
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X (cot4/)j,_i(fcr)-y,_i(fcr)), (All) 



X (cot4/^j,+i(fcr)-y,+i(fcr)). (A12) 



The renormalization of a singular potential in coupled 
channels depends on the number of attractive or repul- 
sive eigenchannels at short distances [1^ [13] ■ Here we 
will only consider in detail the kind of regularization 
which appears when the LO interaction is the OPE po- 
tential |22l |. In such a case, we only need to fix one pa- 
rameter to renormalize the phase shifts. This parameter 
is usually selected to be the scattering length associated 
with the a phase shift, d^j, but in principle there is no 
impediment to fix any other low energy observable. For 
fixing a scattering length we take into account that the 
wave function of the zero energy a and /? states behave 
asymptotically as 



(2j-3)!! 1 
rJ-i (2j - 1)!! oa, 
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(A13) 
(A14) 

(A15) 
(A16) 



where Uaj, Cj and apj are the scattering lengths related 
to the 5aj, tj and djSj phases (i.e. 5aj — >■ —aajk^-'~^ , 
Cj CjP, SjSj -ajSjk'^^'^^), see Ref. for details. 
As previously suggested, we fix Uaj to the desired value, 
leaving Cj momentarily free, and integrate the zero en- 
ergy a state downwards from ?' — !■ oo to r = Tc. At the 
cut-off radius we will impose boundary conditions which 
will determine the initial integration conditions for the 
finite energy a and /3 scattering states. 

There are several possible regularization conditions for 
the coupled channels In the present work, we will 

employ the following set 
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where a, p = (5 and with fj a parameter which is set 
to fj = ^y^j+i for j even and fj = .^/iti for j odd. 

The advantage of these particular regularization con- 
ditions is that they facilitate the energy independent nor- 
malization of the wave functions at short distances. We 
define the new normalization as follows 
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1, we obtain 
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As can be seen, the regularization condition given by 
Eq. (IA17[) assures that the normalization factor A^^j is 
well defined independently of whether we use the u or w 
wave function for computing it. Other regularization con- 
ditions do not guarantee this property, thus generating a 
spurious energy dependence of one of the components of 
the wave function at finite cut-off radii. 
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